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ABSTRACT 

Observations with the Hinode space observatory led to the discovery of predominantly horizontal magnetic 
fields in the photosphere of the quiet internetwork region. Here we investigate realistic numerical simulations 
of the surface layers of the Sun with respect to horizontal magnetic fields and compute the corresponding 
polarimetric response in the Fe I 630 nm line pair. We find a local maximum in the mean strength of the 
horizontal field component at a height of around 500 km in the photosphere, where it surpasses the vertical 
component by a factor of 2.0 or 5.6, depending on the initial and boundary conditions. From the synthesized 
Stokes profiles we derive a mean horizontal field component that is, respectively, L6 and 4.3 times stronger 
than the vertical component. This is a consequence of both the intrinsically stronger flux density of, and the 
larger area occupied by the horizontal fields. We find that convective overshooting expels horizontal fields to 
the upper photosphere, making the Poynting flux positive in the photosphere, while this quantity is negative in 
the convectively unstable layer below it. 

Subject headings: Sun: photosphere — Sun: magnetic fields — MHD — polarization — turbulence 



1. EVITRODUCTION 

Recent observations with the spectropolarimeter of the So- 
lar Optical Telescop e (SOT) o nboard the Hinode space ob- 
servatory (Kosugi et al. 2007) indicate that the quiet inter- 
network region (the inner regions of supergranular cells of 
the quiet Sun) harbors a photospheric magnetic field whose 
mean flux density of the horizontal component considerably 



surpasses that of the vertical component (Lites et al. 2007 
Orozco Suarez et aL]|2007[ Lites et al. 20081. According to 



these papers, the vertical fields are concentrated in the inter- 
granular lanes, whereas the stronger, horizontal fields occur 
most commonly at the edges of the bright granules, aside from 
the vertical fields. In a gravitationally stratified atmosphere, 
vertical magnetic flux concentrations naturally develop a hor- 
izontal component as they expand with height in a funnel- 
like manner Indeed, Rezaei et al. ( 2007 ) found funnel shaped 



magnetic elements in the internetwork from the same Hinode 
data. However, the newly discovered horizontal fields also oc- 
cur apart from vertical flux concentrations and seem to cover 
a larger surface fraction than the vertical fields. 

Regarding numerical simulations, |Grossmann-Doerth et alTj 
([i998 ), note: "we find in all simulations also strong horizon- 
tal fields above convective upflows", and Schaffenberger et al. 
( |2005[ |2006[ ) find frequent horizontal fields in their three- 
dimensional simulations, which they describe as "small-scale 
canopies". Also the 3-D simulations of Abbett (2007 1 display 
"horizontally directed ribbons of magnetic flux that perme- 
ate the model chromosphere", not unlike the figures shown 
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by Schaffenberger et al.| P006| l. More recently, |Schussler| 
& Vogler '{ 2008| l find in a three-dimensional surface-dynamo 
simulation "a clear dominance of the horizontal field in the 
height range where the spectral lines used for the Hinode ob- 
servations are formed". 

Here we report on the analysis of existing and new three- 
dimensional magnetohydrodynamic computer simulations of 
the internetwork magnetic field aiming at following questions: 
Does a realistic simulation of the surface layers of the Sun in- 
trinsically produce horizontal magnetic fields and can their 
mean flux density indeed surpass the mean flux density of the 
vertical field component? What is the polarimetric signal of 
this field and how does it compare to measurements with Hin- 
ode? 

In the following we explain in Sect. 2 the details of two 
simulations and present results that answer the first two of the 
above questions in Sect. 3. In Sect. 4, we synthesize Stokes 
profiles and compare them to measurements from Hinode. 
Conclusions follow in Sect. 5. 

2. TWO SIMULATION RUNS 

We have carried out two runs, run vlO and run h20, which 
significantly differ in their initial and boundary conditions 
for the magnetic field. This enables us to judge the robust- 
ness of our results with respect to magnetic boundary con- 
ditions. Both runs are carried out within a common three- 
dimensional computational domain extending from 1400 km 
below the mean surface of optical depth Tc = 1 to 1400 km 
above it. With this choice we ensure that the top boundary 
is located sufficiently high for not to unduly tamper the at- 
mospheric layers that are in the focus of the present investi- 
gation, in particular the formation layers of the spectral lines 
used in polarimetric measurements with Hinode. The hori- 
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Fig. 1 . — Horizontal (solid) and vertical (dashed) absolute field strength as 
functions of height for run h20 (heavy curves) and for run vlO (light curves). 



zontal dimensions are 4 800 km x 4 800 km, corresponding to 
6.6" X 6.6" on the solar disk. With 120'^ grid cells, the spa- 
tial resolution in the horizontal direction is 40 km, while in 
the vertical direction it is 20 km throughout the photosphere 
and chromosphere. Both runs have periodic lateral bound- 
ary conditions, whereas the bottom boundary is open in the 
sense that the fluid can freely flow in and out of the computa- 
tional domain subject to vanishing total mass flux. The upper 
boundary is "closed", i.e., a reflective boundary is applied to 
the velocity. 

Run V 10 starts with a homogeneous, vertical, unipolar mag- 
netic field of a strength of 1 mT superposed on a previously 
computed, relaxed model of thermal convection. After re- 
laxation, fields of mixed polarity occur throughout the pho- 
tosphere with an area imbalance of typically 3:1 for fields 
stronger than 1 mT. The magnetic field in run vlO is con- 
strained to have vanishing horizontal components at the top 
and bottom boundary but lines of force can freely move in the 
horizontal direction. Although this condition is quite stringent 
for the magnetic field near the top boundary, it still allows the 
field to freely expand with height through the photospheric 
layers. The mean vertical net magnetic flux density remains 
1 mT throughout the simulation. These inital and boundary 
conditions might actually be more appropriate for the simula- 
tion of network magnetic fields because of preference for one 
polarity and the vertical direction. 

Run h20 starts without a magnetic field but upwellings that 
enter the simulation domain across the bottom boundary area 
carry horizontal magnetic field of a uniform strength of 2 mT 
and of uniform direction parallel to the x-axis with them. Out- 
flowing material carries whatever magnetic field it happens to 
have. Thes e boundary conditions are the same as used by 
Stein & Nordlund, ( ,2006^ . They are appropriate when flux as- 
cends from deeper layers of the convection zone, carried by 
convective upflows. Starting from a relaxed model of thermal 
convection, magnetic field steadily spreads into the convective 
layer of the simulation domain and after 600 s slowly begins 
to expand throughout the photosphere, growing in mean ab- 
solute strength. Reflective conditions apply to the field at the 
top boundary, resulting in dB^^y/dz ~ 0, = 0. 

The magnetic energy in the box steadily increases because 
convective plasma motion strengthens the magnetic field. Af- 
ter a time of about 2.45 h an equilibrium value in magnetic 
energy seems to establishes itself when the mean absolute ver- 
tical field strength near the surface of optical depth Tc = 1 is 
approximately 1 mT. Events of convective plumes pump mag- 
netic fields in the downward direction out of the domain so 
that the mean Poynting flux at the lower boundary is negative, 
pointing out of the box. 



Runs vlO and h20 have been carried out with an extended 
version of the computer code CO^BOLD^ that includes mag- 
netic fields. The code solves the coupled system of the equa- 
tions of compressible ideal magnetohydrodynamics in an ex- 
ternal gravity field taking non-local radiative transfer into ac- 
count. For the present runs, frequency-independent opacities 
are used, which are also used for computing the continuum 
optical depth t^. The multidimensional problem is reduced 
to a sequence of 1-D sweeps by dimensional splitting. Each 
of these 1-D problems is solved with a Godunov-type finite- 
volume scheme using an approximate Riemann solver modi- 
fied for a realistic equat ion of state and gravity. Details of the 
method can be found in Schaffenberger et al. ( 2005 2006 1. 



3 . STRUCTURE AND DEVELOPMENT OF THE HORIZONTAL 
MAGNETIC FIELD 

Figure [T] shows the horizontally and temporally averaged 
absolute vertical and horizontal magnetic field strength as 
functions of height for both runs. In run h20, the mean hori- 
zontal field strength, ( \Jb'^ + B'^), is larger than the mean 

strength of the vertical component, throughout the 

photosphere and the lower chromosphere: in run vlO this is 
the case in the height range between 250 km and 850 km. 
It shows a local maximum close to the classical temperature 
minimum at a height of around 500 km, where it is 5.6 times 
stronger than the mean vertical field in case of run h20. The 
horizontal fields also dominate in the upper photosphere of 
run vlO for which case one might expect the initial state and 
boundary condition to favor the development of vertical fields 
rather than horizontal ones. There, the ratio (i?hor) / (^vcr) at 
the location of maximum (i?hor) is 2.5. 

For the second half of the h20 time series and in a horizontal 
section at a height of mean optical depth — 1, 14.2 % of the 
total area is covered by horizontal fields stronger than 5 mT, 
while this fraction is 5.1 % for the vertical fields surpassing 
5 mT. At the height of 200 km in the photosphere and a thresh- 
old of 2 mT the average area fractions are 25.8 % and 6.2 %, 
respectively. Thus, fields with a horizontal component larger 
than a given limit in strength occupy a significantly larger sur- 
face area than fields with a vertical component exceeding this 
limit. This is a second reason (after inherent strength) why the 
measured mean flux density of the horizontal field component 
may exceed that of the vertical one. 

Figure |2] (left) shows for a typical time instant in the sec- 
ond half of run h20 the horizontal field strength (colors) on 
the surface of continuum optical depth Tc = 0.3. Superim- 
posed on the colors are contours of 2 mT of the vertical field 
strength, where solid and dashed contours indicate opposite 
polarity. To the right and the lower right of the image center, 
{x, y) = (3.1 , 2.5) and (3.7 , 1.6), as well as in the middle 
close to the front side, (2.4 , 0.3), we can see a frequently oc- 
curring event consisting of a "ring" of horizontal field. It starts 
to appear as a patch filled with horizontal field like to the left 
front side, (a;, y) — (0.8 , 0.3), subsequently expanding to 
become a ring. The ring can also be seen in the vertical field 
component, where opposite halves of it have opposite polar- 
ity as is visible from the indicated 2 mT contours. This pat- 
tern arises from horizontal magnetic field that is transported 
to the surface by vigorous upflows. The field is anchored in 
the downdrafts at the edges of a granule or in the weaker up- 
wellings of a granule interior. Here, the field is most concen- 
trated and hence not only the vertical but also the horizontal 

2 www.astro.uu. serbf/co5bold_main. html 



The horizontal internetwork magnetic field 
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Fig. 2. — Left: Horizontal field strength on the surface of continuum optical depth Tc = 0.3. The black curves refer to contours of 2mT vertical field strength, 
where solid and dashed contours have opposite polarity. Right: Map of the continuum intensity a 630 nm. 



field is strongest there. As can be seen when comparing to 
the continuum intensity image. Figure [2] (right), a ring does 
often not enclose a full granule but only a part of it so that 
part of the vertical magnetic field occurs within the granule. 
The horizontal field between the crescents of vertical fields of 
opposite polarity covers part of the granule like a cap form- 
ing a small-scale canopy ( Schaffenberger et al. 2005 2006). 



Events of r elated top ograpy were observed by ,Centeno et al.| 
( |2007l ) and |Ishikawa et al. (2008 ). 

As this horizontal field is pushed into the stable layers of the 
upper photosphere by the overshooting convection, it stops 
rising for lack of buoyancy, neither are there vigorous down- 
flows that would pump it back down again. Hence, convective 
flow and its overshooting act to expel magnetic flux from the 
granule interior to its boundaries, i.e., not only to the inter- 
granular lanes but also to the upper layers of the photosphere. 

The surface of Tc ~ 1, which separates the convective 
regime from the subadiabatically stratified photosphere, also 
acts as a separatrix for the verticaUy directed Poynting flux, 
Sz, where 

S=^{Bx{vxB)). (1) 

This can be seen in Fig. [3] (top), which displays the horizon- 
tally averaged Sz as a function of height in the atmosphere 
and time for the second half of run h20. The conspicuous dark 
streaks in the lower part of the diagram mark events of down- 
flow plumes that carry horizontal magnetic field with them, 
giving rise to (Sz) < 0. Differently in the photosphere, where 
(Sz) stays mainly positive (bright), due to the transport of hor- 
izontal fields in the upward direction. They become deposited 
in and give rise to the distinct layer of enhanced horizontal 
fields in the upper photosphere, clearly visible in the middle 
panel of Fig. [3] Here again, both the transport of horizontal 
fields downwards in the convection zone and upwards in the 
photosphere are visible. The bottom panel of Fig.|3]shows the 
mean vertical field strength that monotonically decreases with 
height at all times. 

4. COMPARISON WITH RESULTS FROM THE HINODE SPACE 
OBSERVATORY 

For a reality check we compare the Zeeman measure- 
ments from the Hinode spectropolarimeter with the synthe- 
sized Stokes profiles of both 630 nm Fe I spectral lines of the 
two simulation runs. Profiles were computed with the radia- 
tive transfer code SIR (Ruiz Cobo & del Toro Iniesta 1992[ 
Bellot Rubio|2003" i along vertical lines of sight (disk center) 
with a spectral sampling of 2 pm. We then applied a point 



spread function (PSF) to these 'virtual observations': the the- 
oretical, diffraction-limited PSF of SOT as well as two other 
non-ideal PSFs that take addition al stray-light into account, 
all evaluated at A = 630 nm (see Wedemeyer-Bohm (2008]) 
for details). The following results refer to a PSF obtained 
by convolution of the ideal PSF with a Voigt function with 
7 = 5.7" X 10^"^ and cr = 8" x lO^'^, derived from eclipse 
data. 



For a faithful comparison with the results of Lites et al. 
( 2008 ), we subject the synthetic profiles to the same procedure 
for conversion to apparent flux density^ as was done by these 
authors. Thus, we obtain calibration curves for the conversion 
from the wavelength integrated polarization signals Vtot and 
Qtot to the apparent longitudinal and transversal magnetic 
flux densities |-Bappl ^i^d -^app' respectively. Equally, Qtot is 
the resulting (3-profile after transformation to the "preferred- 
frame azimuth" in which the -i-(5-direction is parallel to the 
projection of the magnetic field vector on the plane of sky, 
when C/ w 0. 

Having the calibration curves, we derive spatial and tem- 
poral averages for the transversal and longitudinal appar- 
ent magnetic flux densities, -Bjpp and |-Bappl of respectively 
21.5 Mxcm^^ and 5.0 Mxcm^^ for the second half of run 
h20 and 10.4 Mxcm"^ and 6.6 Mxcm~^ for run vlO. Thus, 
the ratio r = (Bj^^) / {\Bl; \) = 4.3 in cas e of run h20 
and 1.6 in case of run vlO. |Lites et aL] ( |2008| ) obtain from 
Hinode SP data (|Sjpp|) = 55 Mxcm"^^ (l^appl) = 
11 Mxcm"^ resulting in r = 5.0. 

While (Bhor)/(Svcr) = 5.6 and 2.5 for run h20 and vlO, 
respectively, at the location of maximum (-Bhor), the above 
quoted lower ratios result because the main contribution to the 
Stokes signals does not come from this height but rather from 
the low photosphere, where the two components differ less 
(see Fig. [T]). At full spatial resolution, i.e., without applica- 
tion of the PSF, we obtain from the syntesized Stokes data of 
runh20B'^ = 24.8 Mxcm^^ and \B^\ ^ 8.8 Mxcm^^, thus 
r = 2.8. The higher ratio r when applying the PSF results be- 
cause of apparent flux cancellation within a finite resolution 
element. The vertical component is more subject to this effect 
than the horizontal one because of its smaller spatial scale and 
higher intermittency. This indicates that the predominance 
of the horizontal component decreases with increasing spatial 
resolution and that spatial resolution is a fundamental param- 
eter to take into account when interpreting measurements of 

' Apparent because finite spatial resolution may mask the true flux density 
through cancellation of opposite polarization. 
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Fig. 3. — Vertically directed Poynting flux, Sz, horizontal magnetic flux density, (i?h). ™d vertical absolute magnetic flux density, (|i3z|) as functions of 
height and time from run h20. All quantities are averages in horizontal planes of the three-dimensional computational box. The temporal average of (Sz) is 
maximal 7.4 X 10^ Wm"^ (^t 200 km) and minimal -5.2 X 10* Wm-^ (at -800 km). 

to form a layer of horizontal field of enhanced flux density, 
reaching up into the lower chromosphere. Below the surface 
of Tc w 0.1, convective plumes pump horizontal magnetic 
field in the downward direction. Hence, this surface acts as 
a separatrix for the Poynting flux, which is mainly directed 
upwards above it and in the downward direction below it. 

The response of this field in linear and circular polariza- 
tion of the two neutral iron lines at 630 nm yields a ratio 
(Bjpp)/(|B^pp|) = 4.3 in case of run h20 (which, accord- 
ing to Sect. 2, we deem better to represent the conditions of 
internetwork regions). This is close to the measurements of 
Lites et al. (2008 ), which indicate a factor of 5. Errors may 
come from the straylight produced by the spectrograph and 
polarization optics that was not taken into account with our 
PSF, the difference in mean absolute flux density (run h20, 
has only about half the measured value [see Sect. |4]), the 
frequency-independent treatment of radiative transfer, lacking 
spatial resolution, but also natural fluctuations. The predom- 
inance of the horizontal component, may possibly only exist 
on a scale comparable to or less than the spatial resolution of 
SOT. At full spatial resolution of the simulation we obtain a 
ratio of 2.8 instead of 4.3. 



field inclinations 

density for the fie d inclination at full spatial resolution of the 
simulation shows on the surface Tc — 0.01 a flat (isotropic) 
distribution in the range ±50° from the horizontal direction. 

5. CONCLUSIONS 

We have carried out two simulations of magnetoconvection 
in the surface layers of the quiet internetwork region of the 
solar atmosphere. The simulations greatly differ in their ini- 
tial state and boundary conditions for the magnetic field, but 
otherwise they both equal faithfully reproduce properties of 
normal granulation (as the magnetic field is weak). The top 
boundary is placed in the middle chromosphere (at a height of 
1400 km) far away from the photospheric layers. 

Both simulations intrinsically produce a horizontal mag- 
netic field throughout the photosphere and lower chromo- 
sphere with a mean field strength that exceeds the mean 
strength of the vertical field component at the same height 
by up to a factor of 5.6. The strength of the horizontal field 
component shows a local maximum close to the classical tem- 
perature minimum near 500 km height (which largely escapes 
measurements with the Fe I 630 nm line pair). Fields with 
a horizontal component exceeding a certain limit in strength 
occupy a significantly larger surface area than fields with a 
vertical component exceeding this limit. 

This horizontal field can be considered a consequence of the 
flux expulsion process ( Galloway & Weiss|1981[ l: in the same 
way as magnetic flux is expelled from the granular interior to 
the intergranular lanes, it also gets pushed to the middle and 
upper photosphere by overshooting convection, where it tends 
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